Generic absorbing transition in coevolution dynamics. 
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We study a coevolution voter model on a complex network. A mean-field approximation reveals 
an absorbing transition from an active to a frozen phase at a critical value p c = that only 
depends on the average degree p of the network. In finite-size systems, the active and frozen phases 
correspond to a connected and a fragmented network respectively. The transition can be seen as 
the sudden change in the trajectory of an equivalent random walk at the critical point, resulting in 
an approach to the final frozen state whose time scale diverges as r ~ \p c — p| _1 near p c . 
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The dynamics of collective phenomena in a system of 
interacting units depends on both, the topology of the 
network of interactions and the interaction rule among 
connected units. The effects of these two ingredients on 
the emergent phenomena in a fixed network have been ex- 
tensively studied. However, in many instances, both the 
structure of the network and the dynamical processes on 
it evolve in a coupled manner [H, Hj]. I n particular, in 
the dynamics of social systems (Refs. 0, H,H, [B| and ref- 
erences therein), the network of interactions is not an 
exogenous structure, but it evolves and adapts driven by 
the changes in the state of the nodes that form the net- 
work. In recent models implemen ting this type of coevo- 
lution dynamics @, 0, % H, 0, @, H, M, EL El a transition 



is often observed from a phase where all nodes are in the 
same state forming a single connected network to a phase 
where the network is fragmented into disconnected com- 
ponents, each composed by nodes in the same state [l3j . 

In this paper we address the question of how generic is 
this type of transition and which is the mechanism behind 
it. For this purpose, we introduce a minimal model of in- 
teracting binary state nodes that incorporates two basic 
features shared by many models displaying a fragmenta- 
tion transition: (i) two or more absorbing states in a fixed 
connected network, and (ii) a rewiring rule that does not 
increase the number of links between nodes in opposite 
state. The state dynamics consists of nodes copying the 
state of a random neighbor, while the network dynamics 
results from nodes dropping their links with opposite- 
state neighbors and creating new links with randomly 
selected same-state nodes. This model can be thought 
as a coevolution version of the voter model F| in which 
agents may select their interacting partners according to 
their states. It has the advantage of being analytically 
tractable and allows a fundamental understanding of the 
network fragmentation, explaining the transition numer- 
ically observed in related models [EL d, 0, [l^, 11, H3 |. 
The mechanism responsible for the transition is the com- 
petition between two internal time scales, happening at 
a critical value that controls the relative ratio of these 
scales. 

We consider a network with N nodes. Initially, each 



node is endowed with a state +1 or —1 with the same 
probability 1/2, and it is randomly connected to exactly 
p neighbors, forming a network called degree- regular ran- 
dom graph. In a single time step (see Fig. [I]), a node % 
with state S% and one of its neighbors j with state Sj are 
chosen at random, then: 

1. if Si = Sj nothing happens. 

2. if Si =/= Sj, then with probability p, i detaches its 
link to j and attaches it to a randomly chosen node 
a such that S a = Si and a is not already connected 
to i; and with probability 1 — p, i adopts j's state. 

The rewiring probability p measures the rate at which the 
network evolves compared to the rate at which the states 
of the nodes change; the extreme values correspond to a 
fixed network (p = 0), and to only rewiring {p = 1). 

Link Dynamics. The evolution of the system can be 
described by the densities of two different types of links: 
links connecting nodes with different states or active links 
and links between nodes in the same state or inert links. 
Note that an update (either rewire or copy) only occurs 
when an active link is chosen. 

In Fig. [1] we describe the possible changes in the global 
density of active links p and their probabilities in a single 
time step, when a node of degree k is chosen. We denote 
by n the number of active links connected to node i be- 
fore the update. With probability n/k an active link i—j 
is randomly selected. Then with probability p the link 
i — j is rewired and becomes inert (link i — a), giving a 
local change of active links An = — 1 and a global density 
change of Ap — — -J^, where pN/2 is the total number 
of links, fi — (k) — J^k^Pkit) is the number of links 
per node or average degree and Pk (t) is the node degree 
distribution at time t. (Even though Pk(t) depends on 
time, given that the network is constantly evolving, p is 
constant because the total number of links is conserved 
at each time step. Furthermore, simulations show that 
Pkit) has a narrow shape with a maximum at pi, as ex- 
pected from the random nature of the rewiring). With 
probability (1 — p) node i flips its state changing the state 
of links around it from active to inert and vice- versa, and 
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FIG. 1: Update events and the associated changes in the 
density of active links p and the link magnetization m — 
p++ — when two neighbors i and j with states Si — s 
and Sj — — s are chosen (s = ±1). 
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FIG. 2: Stationary density of active links p s and the two types 
of inert links p++ and p__ vs the rewiring probability p as 
described by the mean-field theory for a network with average 
degree p — 4. The critical point p c separates an active from 
a frozen phase. 



leading to An = k — 2n and Ap — 2 ( fc ~^") ■ Assembling 
these factors, the change in the average density of active 
links in a single time step of interval dt = -h is described 
by the master equation 
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where H n ,k is the probability that n active links are con- 
nected to a node of degree k, and (n)t and {n 2 )k are the 
first and the second moments of B n< k respectively. In a 
mean-field (MF) spirit, in which the system is considered 
to be homogeneous, we approximate the probability that 
a given link is active by the average density p. Then B n k 
becomes the binomial distribution with first and second 
moments (n)k = pk and (n 2 )k = pk + p 2 k(k — 1). Re- 
placing these expressions in Eq. (fTJ) we obtain a closed 
equation for the time evolution of p 

f = ^[(l-p)(/x-l)(l-2p)-l]. (2) 

Equation ([2|) has two stationary solutions. For p < p c , 
the stable solution is 



£(?) = 



(i- p )(p-i)-i 

2(l-p)( M -l) : 



(3) 



corresponding to an active steady-state with a constant 
fraction of active links in the system, while for p > p c the 
stable solution p s = corresponds to an absorbing state 
where all links are inert. Thus, the MF approach predicts 
an absorbing transition (see Fig. [2]) from an active to a 
frozen phase at a critical value 



p-1 



(4) 



A stationary state in the active phase, characterized by a 
density £ (p) = 2(i-p) °^ ac ti ye links, is composed by links 
continuously being rewired (evolving network) and nodes 
flipping their states. For p = (original voter model), the 
value £(0) = 2 (^-i) a S rees very well with the numerical 
values of the voter dynamics in different random graphs 
3, ID]. In the frozen phase, final states correspond to 
a fixed network where connected nodes have the same 
state and no more evolution is possible. The transition 
is continuous, with the order parameter p s changing con- 
tinuously at p c . Close and below the transition point, p B 
scales as (p c —p), thus the MF critical exponent is 1. 
In order to obtain an insight about the structure of 

the network in both phases we introduce (p ) as 

the density of links connecting two nodes with states +1 

(—1). It can be shown that p ++ (p ) is related to p 

and the density a + (cr_) of + (— ) nodes by 



p ++ = a + - p/2; p__ = cr_ - p/2. 



(5) 



Due to the conservation of the ensemble average of <j + 
and <7_ for the voter model dynamics, we have that 

cr+ = cr_ = 1/2, and therefore, p++ = p = ^(1 — p). 

In the active phase the continuous rewiring of links keeps 
the network connected in a single component, i.e., a set 
of connected nodes. But, in the frozen phase only in- 
ert links are present and in the same proportion (p++ = 
p__ = 1/2), thus we expect the formation of two large 
disconnected components with opposite state (see Fig. [5]). 
Therefore, the MF description reveals a fragmentation 
transition in the stationary structure of the network, as- 
sociated with the absorbing transition at p c . 

Final states in a finite system. The previous MF ap- 
proach predicts a transition in the limit of an infinite 
large network. However, for any value of p, due to fluc- 
tuations, a finite size network eventually reaches an ab- 
sorbing state composed by inert links only. We studied 
the structure of the network in the final state by per- 
forming numerical simulations of the dynamics starting 
with a degree-regular random graph with connectivity 
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FIG. 3: (a) Average relative size of the largest network com- 
ponent S and (b) absolute value of the link magnetization m 
vs p in the final frozen state, (c) Average convergence time r 
per system size N vs p. Inset: scaling of r for p > p c ~ 0.38, 
indicating that r ~ ~ ln(z) with z = p(p — p c )N. The solid 
line has slope —1. (d) Stationary density of active links in 
surviving runs p surv . Inset: average time evolution of p. The 
averages are over 10 4 realizations of networks with p — 4 and 
sizes N = 250 (circles), 1000 (squares) and 4000 (diamonds). 



fi = 4 and letting the system evolve until it was frozen. 
In Fig. Ela) we plot the average size of the largest net- 
work component S in the final configuration for networks 
with N = 250, 1000 and 4000 nodes. We observe that S 
is very close to N for values of p below a transition point 
p c ~ 0.38, indicating that the network forms a single 
component fl7| . Above p c the network gets disconnected 



into two large components and a set of components of 
size much smaller than N, giving a value of S ~ N/2. 

To compare simulations with MF results, we calcu- 
lated the stationary density of active links in surviving 
runs p surv . As we show in Fig. [3Jd), p surv monotonically 
decreases with p, becoming sharper with increasing sys- 
tem and indicating a transition from an active to a frozen 
phase as predicted by the MF theory. In the active phase, 
p surv reaches a steady value larger than zero and inde- 
pendent on the system size N, while in the frozen phase 
/3 surv vanishes in the thermodynamic limit. The critical 
point for the active-frozen transition p c ~ 0.38 calculated 
from Fig. [3jd) is roughly the same as for the fragmen- 
tation transition (Fig. |3£a)), suggesting that the active 
and frozen phases observed in infinite large systems corre- 
spond to the connected and disconnected phases respec- 
tively in finite systems. The MF critical point p c = 2/3 
calculated using Eq. Q with p, = 4 differs from the nu- 
merical value p c ~ 0.38 (Figs.0Ja,d)) due to correlations 
appearing in the rewiring process. These correlations, 
that are not taken into account in the analytical approx- 
imation, make the moments (n)k and (n 2 )^ different from 
the analytical values pk and pk + p 2 k(k — 1) respectively. 
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FIG. 4: Typical trajectories of the random walk for a network 
of size N — 10 4 and p — 4. The upper (p = 0) and lower 
(p = 0.35) parabolas are the trajectories for rewiring rates 
below the transition point p ~ 0.38, while the quasi-vertical 
line is for p — 0.4. Insets: Time evolution of the density of 
active links p (right) and the ratio between p and the product 
<7+<7_ (left) in a single realization for different values of p. 



These deviations have the overall effect of decreasing the 
observed critical point respect to the theoretical one. 

Approach to the absorbing states. So far, we have 
shown that a finite network under the coevolving dy- 
namics experiments a fragmentation transition as the 
rewiring rate is increased. We now unveil the mecha- 
nism of this transition by studying the evolution of the 
system to the frozen state. 

We represent the state of the system as a point (m, p) 
in the 2 dimensional space, where the coordinates are 
the link magnetization m = p++ — p and the den- 
sity of active links respectively. When a node of de- 
gree k connected to n active links is chosen, the possible 
changes in m and p and their respective probabilities are 
those described in Fig. [TJ In the (m, p) space, the sys- 
tem undergoes a random walk (RW) inside the triangle 
< p + \m\ < 1, whose boundaries follow from the con- 
straint relation p__ + p ++ + p = 1. The system reaches 
an absorbing configuration and stops evolving when the 
RW hits either one of the fixed points (—1,0) or (1,0) 
(all nodes in state — or + respectively) or a point on the 
fixed line p — (frozen mixture of — and + nodes). At 

the point (—1, 0) ((1, 0)) only (++) links are present, 

the network is composed by a giant component, and the 
system is in the connected phase. Points on the line p = 
and close to point (0, 0) correspond to a frozen network 

with similar number of H — h and links arranged in two 

large + and — components (disconnected phase). 

In Fig. [4] we plot trajectories of the RW in one re- 
alization for different values of p. For p < p c ~ 0.38, 
the motion of the RW has two stages. In the first and 
very short stage, the RW travels along the m ~ axis 
from the starting point (m ~ 0,p ~ 1/2) to the point 
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(m ~ 0, p ~ £ p ) that corresponds to the steady state 
p s = in infinite large systems (see right inset of FigH}. 
In the second and long stage, the RW diffuses on the m 
direction, corresponding to the fluctuations of p down the 
steady state p s = £ p , until it hits either point m = — 1 
or m = 1 (see Fig. EJb)). We observe in Fig. Q] that 
the motion of the RW is not completely random but its 
trajectory fluctuates around a curve described by 

P P (m) = f P (l ~ m2 )- ( 6 ) 
The origin of this relation is that both p and m can be ex- 
pressed as functions of er + . As we observe in the left inset 

fluctuates 



of Fig. [U during one realization the ratio 
around the constant value 4£ p . Then, using Eq. ([5]) we 

obtain p = 4£ p er_|_(l— er+) and m = p ++ —p = 2<r + — 1, 

from where we arrive to Eq. ([6]) by eliminating cr+. For 
p > p c , the bias to the p = line makes the RW hit a 
point close to the origin (see the p = 0.4 trajectory in 
Fig. |4|). Simulations show that for a fixed value p < p c , 
the amplitude of the fluctuations of the RW's trajectory 
around its mean value p p (m) vanishes as N increases. 
Thus, increasing N has the effect of increasing the prob- 
ability that the RW reaches one of the end points m = ±1 
before it hits the line p = 0, and therefore, most realiza- 
tions end in a single component. Eventually, in the large 
N limit, the RW has three absorbing points: either point 
(—1, 0) or (1, 0) (single component network) whenp < p c , 
and point (0, 0) (two components network) when p > p c . 

Convergence times. A magnitude of interest is the av- 
erage time t to reach an absorbing state. For p < p Cl 
the m-coordinate of the walker performs a Id symmetric 
random walk with an average jumping interval and its 
probability that scale as 1/N and p ~ £ respectively. To 
reach one of the end points m = ±1 the RW needs to at- 
tempt an average of iV 2 /£ steps, and given that the time 
increases by 1/N in each attempt, we find that r ~ N/£. 
From Eq. ([2]) for p > p c , p decays as pit) ~ — £e 4 ^*/ M . 
The system freezes at a time r for which p(r) ~ 1/pN, 
then r ~ — ^| ln(— fi£N). Using the MF approximation 
£(p) ~ [Pc~p) close to pc, we obtain that r ~ iV(p c — p) -1 
as p ^ p~ and r ~ (p - Pc) -1 ln[/i(p - p c )-/V] as p — > p+, 
thus the convergence to the final state slows down at the 
critical point (see FigJ^c)). 

Summary and conclusions. In summary, the coevolu- 
tion mechanism on the voter model induces a fragmenta- 
tion transition that is a consequence of the competition 
between the copying and the rewiring dynamics. In the 
connected active phase, the system falls in a dynami- 
cal steady state with a finite fraction of active links. The 
slow and permanent rewiring of these links keeps the net- 
work evolving and connected until by a finite-size fluctu- 
ation the system reaches the fully ordered state (all nodes 
in the same state) and freezes in a single component. In 
the frozen phase, the fast rewiring dynamics quickly leads 
to the fragmentation of the network into two components, 
before the system becomes fully ordered. 



The similarity between the mean-field equation for the 
density of active links in the coevolution voter model 
(Eq. [2|) and the one for the density of infected sites in the 
contact process [3], suggests that our model could be- 
long to the Directed Percolation universality class. How- 
ever, both models seem not to be equivalent given that 
our model possesses many absorbing states (any point 
on the p = line of Fig. HJ , unlike the contact process 
where there is a single absorbing state characterized by 
the absence of infected sites. We believe that our results 
provide a new insight in the ongoing discussion about 
models with infinitely many absorbing states [20| . 
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